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Abstract 

Harmonic sums and their generalizations are extremely useful in the evaluation of 
higher-order perturbative corrections in quantum field theory. Of particular interest 
have been the so-called nested sums, where the harmonic sums and their generaliza- 
tions appear as building blocks, originating for example from the expansion of gener- 
alized hypergeometric functions around integer values of the parameters. In this Letter 
we discuss the implementation of several algorithms to solve these sums by algebraic 
means, using the computer algebra system Form. 
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Version: 1.0 
Catalogue number: 
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1. Introduction 



Symbolic summation amounts to finding a closed-form expression for a given sum or 
series. Systematic studies have been pioneered by Euler [1], and for specific sums, ex- 
act formulae have been known for a long time, series representations of transcendental 
functions being a prominent example. Today, general classes of sums, for example 
so-called harmonic sums, have been investigated (see e.g. Refs. [2]) and symbolic 
summation has further advanced through the development of algorithms suitable for 
computer algebra systems. Here, the possibility to obtain exact solutions by means of 
recursive methods has lead to significant progress, for instance in the summation of 
rational or hypergeometric series, see e.g. Ref. [3]. 

In quantum field theory, higher-order corrections in perturbation theory require the 
evaluation of so-called Feynman diagrams, which describe real and virtual particles 
in a given scattering process. In mathematical terms, Feynman diagrams are given as 
integrals over the loop momenta of the associated particle propagators. These inte- 
grals may depend on multiple scales and are usually divergent, thus requiring some 
regularization. The standard choice is dimensional regularization, i.e. an analytical 
continuation of the dimensions of space-time from 4 to D, which keeps underlying 
gauge symmetries manifestly invariant. Analytical expressions for Feynman integrals 
in D dimensions may lead to transcendental or generalized hypergeometric functions, 
which have a series representation through nested sums with symbolic arguments. The 
main computational task is then to obtain the Laurent series upon expansion of the 
relevant functions in the small parameter £ = (D — 4)/2. 

It is the aim of the present Letter to discuss the implementation of several algo- 
rithms [4] for these tasks in the computer algebra system Form [5, 6]. The result- 
ing package XSummer has already been used in full-fledged calculations in particle 
physics, for instance in calculating higher-order perturbative corrections in Quantum 
Chromodynamics, see e.g. Ref. [7]. We hope that it may also be useful for a larger 
community, as it exceeds current built-in capabilities of commercial computer algebra 
systems such as Maple or Math EM ATI CA in the expansion of (generalized) hypergeo- 
metric functions. 

To give a concrete example for the kind of problems we aim at, consider the hyperge- 
ometric function which has a series representation for arguments |jc| < 1: 



2F1 (ae, be, l-ce,x)^Y, \ 



= \^-ab\l2{x)£^ + ab{cLi2,{x) + {a + b + c)Si^2{x)]£^ + 0{z'), (1.1) 

where the {az)^ = r(7-|-fl8)/r(a£) are so-called rising factorials (in the literature also 
known as Pochhammer symbols). Here, we have expressed the coefficients of the 
Laurent series in £ through standard polylogarithms Li„ and Nielsen functions S„,p, 
see e.g. Ref. [8]. 
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Our choice of Form for the implementation is based on two reasons, one being that 
Form is extremely fast and flexible when dealing with large expressions. Form allows 
for a very compact notation and is equipped with a pattern matcher well suited to solve 
our problem. The ability to handle large-size expressions (of the order of the computer 
memory) is crucial, as they generally occur at intermediate stages in the quantum field 
theory calculations mentioned above, for instance when the Laurent expansions in 8 
have to be done to very high order. 

The other main motivation for choosing Form is due to the existing Summer pack- 
age [9]. The Summer package in Form is capable of solving nested sums in terms of 
harmonic sums, a feature that has also been used extensively in recent cutting-edge cal- 
culations of structure functions to three loops in Quantum Chromodynamics [10-12]. 
Here, the XSummer package, being capable of handling (multiple) scales, e.g. in 
Eq. (1.1), provides the obvious extension. In the development of XSummer, we have 
also benefited from the fact that some parts of the underlying algorithmic structures 
could be literally taken from Summer. 

A major disadvantage of using Form for XSummer is certainly the lack of internal 
algorithms for particular operations on polynomials, such as factorization. We will 
comment on that in the text. Also, we note that an implementation of the algorithms 
of Ref. [4] within the GiNaC framework [13] is available [14] . 
The outline of the Letter is as follows. In Section 2, we briefly give the basics of 
generalized sums and recall the algorithms of Ref. [4]. In Section 3, we present the 
XSummer package and discuss details of the implementation. Section 4 features an 
extensive set of tests with various sample calculations, including e.g. Ref. [7]. We 
conclude in Section 5. 

2. Harmonic sums and their generalizations 

The basic recursive definition of iS-sums is given by [4] 



These are the basic objects that we will manipulate in the following. Generally, we 
have all \xi\ < 1 in Eq. (2.1). The sum of all nij is called the weight of the sum, while 
the index k denotes the depth. This definition actually includes as special cases the se- 
ries representations of classical polylogarithms, Nielsen functions, as well as multiple 
and harmonic polylogarithms [15-17]. For all Xi = 1, the above definition reduces to 
harmonic sums [1,9, 18] and, if additionally the upper summation boundary n — > oo, 
one recovers the (multiple) zeta values associated to Riemann's zeta-function [2]. 




S{n;mi,...,mk;xi,...,Xk) 



(2.1) 
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An equivalent representation of S-sums reads 



S(n;mi,...,mk;xi,...,Xk) = Y (2.2) 



We note that S-sums are closely related to so-called Z-sums, the difference being the 
upper summation boundary for the nested sums: (4 ~ 1) for Z-sums, ij^ for S-sums, 
see Ref. [4]. One can algebraically convert Z-sums to 5-sums and vice versa [19,20]. 
We rely entirely on 5-sums in our discussions, but nevertheless provide the procedure 
ConvStoZ since, in some cases, Z-sums may be more favourable. 
The iS-sums obey the well-known algebra of multiplication, a straightforward general- 
ization of the results on the multiplication of harmonic sums [9]. The basic formula 
reads 

S{n;mi, ...,mk;xi, ...,Xk) x S{n;m'i, ...,mi;x'i,...,x'i) 

n 

-:M^S{h;m2, ...,mk;x2, ...,Xk)S{ii;m[, ...,m'i;x[,...,x'i) 



h=i h 



n ^12 



+ L -j-rS{i2-,mi, ...,mk;xi, ...,Xk)S{i2\m2, ...,mf,X2, ...,xi) 



.m, 
12=1 «2 

~£ S{i\m2,...,mk;x2,...,Xk)S{i;m2,...,m'i;x'2,...,Xi), (2.3) 

1=1 * 

which works recursively in the depth of the individual sums. The algorithm allows 
the expression of any product of nested sums as a sum of single nested sums, hence 
in a canonical form, which is an important feature for practical applications. The 
underlying algebraic structure in Eq. (2.3) is a Hopf algebra, being realized as a quasi- 
shuffle algebra here, see e.g. Refs. [2, 4] . The algorithm can be implemented very 
efficiently on a computer, see the procedure BasisS in Section 3. 
In our applications, such as in Eq. (1.1), we encounter Gamma-functions, which we 
have to expand in the small parameter 8 before the actual manipulation of the nested 
sums. This proceeds according to the well-known formula for the expansion of the 
Gamma-function around positive integer values. 

Similarly, the expansion of the Gamma-function around negative integer values can be 
reduced to the case in Eq. (2.4) with the help of the following relation (e.g. Ref. [21] 
p. 3), 

r(i+8) ^ ' r(n-8)' ^ ^ 
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which yields 

\n-l 



2.1. Algorithms 

For the manipulation of the 5-sums, we classify four types of transcendental sums. 
These types of sums are dealt with in the algorithms A to D given below. All sums in 
these classes can be solved recursively, i.e. they can be expressed in canonical form. 
The underlying algorithms realize a creative telescoping. They either reduce succes- 
sively the depth or the weight of the inner sum, so that eventually the inner nestings 
vanish and the results can be written in the basis of Eq. (2.1). The procedure gen- 
erally relies on algebraic manipulations, such as partial fractioning of denominators, 
shifts of the summation ranges and synchronization of summation boundaries of the 
individual sums. Another crucial ingredient is, of course, the quasi-shuffle algebra of 
multiplication in Eq. (2.3). 

Basic definition (type A): 

Here we consider sums over / involving only S{i;...), of the form, 

S{n;mi,...,mk;xu...,Xk) = — — J— -5'(/-F&;m2, ...,m*;;x2, , (2.7) 

i=l V i^J 

where we assume that a, b are (non-symbolic) integers. The upper summation limit is 
allowed to be infinity. 



Convolution (type B): 

Here we consider convolutions, i.e. sums over / involving both S{i; ...) and S{n — i; ...), 
of the form, 

n—l. -^i 

[Tr-\^S{i + b\m2,...,m],\X2,...,Xk) 
(=1 ^yi'T'O') 

^ / ^^^^ ,.^f S{n-i + b';m2,...,m'i;x'2,...,Xi)) , (2.8) 

where all a, a', b and b' must be (non-symbolic) integers. Note that the upper summa- 
tion limit is (n — 1) and thus consistent with the defining range of the S-sums. 
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Conjugation (type C): 

Here we consider conjugations, i.e. sums over / involving 

{-iys{i;mi,...;xi,...) 



and a binomial 



of the form, 

( / ) {-iyjr^^^S{i + b-m2,...,mk;x2,...,Xk), (2.9) 



where a^b are (non-symboUc) integers. The upper summation limit should not be 
infinity. Again, the upper summation limit n is consistent with the defining range of 
the binomial. Sums of this type cannot be reduced to S-sums with upper summation 
limit n alone. However, they can be reduced to S-sums with upper summation limit n 
and multiple polylogarithms (which are S'-sums to infinity). 

Binomial convolution (type D): 

Here we consider binomial convolutions, i.e. sums over / involving (— l)'5'(i; ...), 
S{n — i;...) and a binomial, of the form, 

/ fl \ x' 

^T~^-T^~;zzr^i^-^ + b''^^2i---imf,X2,...,Xi). (2.10) 

Here, all a, a', b and b' must be (non-symbolic) integers. Yet again, the upper summa- 
tion limit {n — I) reflects the defining range of the binomial and the S'-sums. As for 
sums of type C, we cannot relate them to S-sums with upper summation limit (n—l) 
alone, but we can reduce them to 5-sums with upper summation limit (n—l) and 
multiple polylogarithms (which are i'-sums to infinity). 



3. The XSUMMER package 

In this section we give a short description of the XSummer package. In particular we 
explain the notation that we use in the package and the main routines that act as a 
front-end to a bunch of smaller routines used internally. At the end of this section we 
also comment briefly on the internal routines, although the user might not want to call 
them directly. 
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3.1. Basic syntax — form follows function 

The XSUMMER package is written using the computer algebra system Form. Unlike 
programs such as MAPLE or Mathematica the program Form provides only a very 
limited set of built-in capabilities. Form is mainly a highly efficient pattern matcher. 
In choosing a syntax/notation for the XSummer package it is therefore important to 
ensure that all basic algorithms can be implemented as simple pattern matching and, 
at the same time, make this pattern matching as simple as possible. 
Generically, we use the function R with an arbitrary number of arguments to denote 
a list of integer parameters. Similarly, X is used to denote a list of symbolic argu- 
ments. Due to the internal limitations of Form with respect to polynomial algebra, it 
is of some advantage to define, for various functions, a distinct multiplicative inverse, 
which helps in bringing expressions to a normal form. Examples are the pairs Gamma 
and InvGamma. Also note, that the summation symbol simply appears as a function 
multiplying all terms that should be summed over. For instance, a sum of the form 
shown in Eq. (2.9), 



would be written as: 

sum(jl,l,n) * bino(n,jl) * pow(-xl,jl) * den(jl+3)^2 * 
S (R(4,7, 1,1) ,X(x2,x3,x4,x5) , jl+2) ; 

Table 3.1 shows the complete set of keywords for the XSummer package. All the 
objects in Table 3.1 are defined in the file declvars.h which should be included when 
using the XSummer package. 

Note that the pow function is reserved for symbolic variables to some power of a 
summation index. Writing 

pow( jl+3,-2) 

instead of 

den{ jl+3) "2 

would not work. This is just due to the way the pattern matching is realized in the 

package. 

Furthermore, it is assumed that the summation variables are always j\,j2, — where the 
outermost sum runs over ji, the next sum runs over 72 and so on. The innermost sum is 
the sum over the highest ji. Upon summation of a nested sum the innermost sum must 
be done first and we then work through to the outer sums. This is essentially done by 
calling the procedure DoSum, which takes as arguments the summation indexes of the 
innermost and the outermost sum. For example: 




(3.1) 



7 



Name Description Example 

bino binomial coefficient bino(n,i) ^ (") 

delta Kronecker delta delta(x) ^ | o 

deltap inverse Kronecker delta deltap(x) ^ ( 1 - delta(x) ) 

den denominator function den(x) — ^ 

ep expansion parameter £ ep — > 8 

epow powers of epsilon epow(n) 8" 

fac factorial function fac(n) — > n\ 

inf parameter for oo inf — > oo 

invfac inverse factorial function invfac(n) ^ 

num numerator function num(x) x 

pow power function pow(x,a) 

sign sign function sign(n) ( 1)" 

sum summation symbol sum(j,i1,i2) ^ ^^JL,^ 

theta theta function theta(x) | o l<o 

z2, z3,. . . values of the zeta-function z2-^ C(2), ■ • • 

Gamma Gamma-function Gamma(x) r{x) 

InvGamma inverse Gamma-function InvGamma(x) tJ-t 

r{A-) 

S(R(mi , . . .),X(^i , . . .),n) 5-sum S(R(mi , . . . ,),X(xi , .),") 

S{n;mi, . . . ,mk;xi, . . . ,Xk) 
Z(R(mi , . . .),X{xi , . . .),n) Z-sum Z(R(mi , . . . ,),X{xi , . . .),n) 

Z{n;mi, . . . ,mk;xi, . . . ,Xk) 

Table 1: Basic objects appearing in the input/output of the XSummer package 
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#call DoSum(3,l) 

would evaluate the sums containing 73, 72. ji (in this order). Note that, when doing the 
sum over a specific ji, all the objects relevant for this sum are dressed internally with 
an additional index / as part of the names of the symbols. For example the sum shown 
in Eq. (3.1) would be converted internally to 

suml(jl,l,n) * binol(n, jl) * powl(-xl,jl) * denl(jl+3)^2 * 
SI (R(4,7,l, 1) ,X(x2,x3,x4,x5) , jl+2) ; 

This is just to simplify the pattern matching. However, in the final result the dressed 
objects should never occur. The following Form script solves a much easier version of 
the example shown above (otherwise the result would be to lengthy to be reproduced 
here): 

#define MAXSUM "1" 
♦define MAXWEIGHT "20" 
♦include declvars.h 
nwrite stat; 

L demo = sum(jl,l,n) * bino(n,jl) * pow{-xl,jl) * den{jl+l) 
* S(R{l,l),X{x2,x3), jl + 1); 

id bino(xl?,x2?) = fac (xl) *invfac (x2) *invf ac (xl-x2) ; 

#call DoSum(l,l) 

print; 

.end; 

The result obtained from running FORM is given by: 

demo = 

+ acc (-1) *pow (1 - xl,l + n)*den( - xl)*den(l + n) 
*S(R(l,l),X(den(l - xl) - den(l - xl) *xl*x2, den (1 - xl*x2)) 

,1 + n)*theta( - 1 + n)+ acc(l)*pow(l - xl,l + n)*den( - xl) 
*den(l + n) *S(R(1, 1) ,X(den(l - xl)- den(l - xl)*xl*x2, 
den{l - xl*x2) - den(l - xl*x2) *xl*x2*x3) , 1 + n)*theta{ - 1 + n) 
+ acc(l)*den( - xl)*theta( - 1 + n)*xl*x2*x3 

} 

The acc function is defined internally with the Form declaration PolyFun to collect 
similar objects together, in particular to accumulate powers of the expansion parameter 
E. Its use is actually described in the FORM manual and the expanded result is simply 
obtained with 
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id acc (x?) 



= x; 



As mentioned earlier the integer parameters of the 5- sums are collected in R, while the 
symbolic arguments are collected in X. Note that only basic simplifications are applied 
to the symbolic arguments. This is due to the fact that Form does not provide built-in 

routines to factorize or normalize expressions. Such procedures must be provided by 
the user and are highly dependent on the problem studied. We will discuss this issue 
in Section 4. 

Our example, converted back to a more human readable notation is given by 
t ( • ')^-^5(7i + l;l,l;x2,X3) = 



1 (l-xi^"+i 



0(n - \){—^^- ""-^ kn+ 1; 1, 1; (1 -xiX2)/(l -xj), 1/(1 -xyxi)) 

-5(n+l;l,l;(l-xiX2)/(l-xi),(l-xiX2X3)/(l-xiX2)) -X2X3|. (3.2) 

Upon replacement n — > the right-hand side of Eq. (3.2) could serve as input to 
another summation, thus nicely illustrating the telescoping character of the recursions 
discussed in Section 2. 

In passing, we note that in some cases finite polynomial sums appear. For example we 
may encounter sums of the types 

n n n—l n—l 

7=1 7=1 7=1 7=1 

Internally they are called xpowsum, powsum, xpowsumi, powsuml. For positive 
integer values of k up to 10 these sums are tabulated in the file declsums.h. Gen- 
erally, for any integer k, these types of sums can easily be obtained with programs 
such as Maple or Mathematica, where the corresponding recursive algorithms are 
implemented, see e.g. Ref. [22]. 

3.2. List of procedures provided byXSuwiwiER 

As briefly mentioned above, we distinguish two sets of procedures: those easily acces- 
sible to the user and those for internal use only. 

3.2.1. Procedures to be called directly 
BasisS 

Express products of S'-sums into single S'-sums of higher weight. 
ConvStoZ 

Convert S-sums into Z-sums. 
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DoSum 




Makej 


|i 


DoTheta 




DoSynch 




SplitDen 


N 


DoTheta 




RegSum 




Sumit 




SimplifyFacGam 




SimplifySumArgs 




SubsS 




BasisS 



DoSynch 




SynchSS 




NormDen 




SplitDen 




Flipj 




Movej 




SynchNumFac 




SynchDenS 




SynchDenFao 




Synch FacS 




SplitDen 




DoTheta 




SynchNumFac 




Flipj 




SynchDenS 




Adjusts urn 



■ SynchSS 
BasisS 



-SumIt 

AlgD 

IAIgC 
AlgBO 
AlgB 
SumPosPow 



■ SynchNumFac 
SubsS 



RegSum "^^^^P 




Flipj 




SynchDenS 




Movej 




Flipj 




Movej 




DoTheta 




SubsS 



AdjustSum 




DoTheta 




SubsS 




Movej 




SubsS 




Movej 




SubsS 



■AlgD 



AlgC 



AlgC 



SimpleArgs 



■ AlgBO 
AlgB 



Figure 3.1: Internal structure of the XSummer package 
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DoSum 

User front-end to the XSummer package. The first parameter denotes the 
index of the innermost summation, the second is the index of the outermost 
sum, which should be summed. For example 

#call DoSum{3,l) 

would do the sums over 73, j2, ji in this order. 

An overview of the internal structure of DoSum is shown in Fig. 3.1. 

3.2.2. Internal procedures 
AdjustSum 

Adjusts the sum boundaries. The argument specifies the index i of the sum to 
be adjusted. 

AlgB 

Implementation of the convolution algorithm (type B) to perform sums of the 
type 

c/- ^ ^ 1'' ■ I II i\ 

—S{i;m2, ...,mk;x2, ...,Xk)- —rS{n - i;m2, ...,m;;x2, , 

for any (integer) values of m/, and mi > 0. The argument specifies the index 
/ of the sum to be done. The special case mi = is handled in AlgBO. 

AlgBO 

Special case of the convolution algorithm (type B), for details see AlgB. 
AlgC 

Implementation of the conjugation algorithm (type C) to perform sums of the 
form 

L " ) (-l)'^'5(«;'W2,...,mA;;x2,...,Xfe) 

for any (integer) values of m; and mi > 0. The argument specifies the index / 
of the sum to be done. The special case mi = is treated in AlgCO. 

AlgCO 

Special case of the conjugation algorithm (type C), for details see AlgC. 
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AlgD 

Implementation of the binomial convolution algorithm (type D) to perform 
sums of the form 



,^1 ( ) ^"^'''/i'^*-'''"^'"'''"'^'^^'"''^*'' 
(n - ip^ 

for any (integer) values of m/,m- > and mi > 0. The argument specifies the 
index / of the sum to be done. The special case mi = is treated in AlgDO. 

AlgDO 

Special case of the binomial convolution algorithm (type D); see AlgD. 
DoSynch 

Procedure to synchronize the arguments of num, den, fac, invfac and S. The 
argument specifies the index i of the sum to be adjusted. Start with synchro- 
nizing products of S functions, then do the combinations 

1. ji and fac, invfac 

2. den and S 

3. den and fac, invfac 

4. S and fac, invfac. 

Finally summation boundaries are synchronized. 
DoTlieta 

Simplifies combinations of theta, delta, deltap and sum functions, defined 
for each sum over jj. The argument specifies the index / of the sum to be 
adjusted (argument / = implies no sum). 

Expand Den 

Expands den function in small parameter ep, i.e. l/{a + be) in terms of 8. 
If the argument provided to Expand Den is then denominators with integer 
values of a are expanded. If the argument is 1, we expand for symbolic a. 

ExpandGam 

Expands Gamma and InvGamma functions in the small parameter ep, i.e. 
r{i + ae) and 1 /r{i + ae) in £, where i and a take integer values. The argument 
specifies the number of the highest sum in the expression. We choose the 
MS-scheme, i.e. exp(— y^ae) = 1, where Je is Euler's constant. 
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Flipj 

Reverses the direction of summation. The argument specifies the index / of 
the sum to be flipped. 

Makej 

Creates the dressed objects, for example the pow function pow (x, j3 ) is con- 
verted to pow3 (x, j3) . The argument specifies the summation index to be 
treated. 

Movej 

Makes a translation of ji + n — > ji. The first argument specifies which index 
/ of the summation variable which should be used. The shift is determined 
by the argument of the function specified by the second argument passed to 
Movej. 

Norm Den 

Converts den functions to normal form, i.e. fixes the sign of the summation 
index appearing in the denominator. The argument specifies the highest sum- 
mation index i to be treated. The hierarchy is such that j\ is more important 
than 72 etc. 

NormalizeGam 

Normalizes products of Gamma and InvGamma functions. 
RegSum 

Performs final regularization of a sum. The argument specifies the index / of 
the sum to be regularized. 

SimpleArgs 

Applies some trivial simplifications. If the argument passed to SimpleArgs is 

these simplifications are applied to the prefactors multiplying functions. If 
the argument passed is C the simplifications are applied to the arguments of S 
and pow functions. 

Simplify 

Calls different simplification and normalization routines to simplify or nor- 
malize the expresssions. 

SimplifyFacGam 

Procedure to simplify products of fac, invfac and Gamma, InvGamma func- 
tions. 

SimplifySumArgs 

Procedure to simplify polynomials in arguments of S functions. This proce- 
dure is user accessible and should be edited if optimization is needed. 
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SplitDen 

Partial fractioning of products of denominators involving the summation vari- 
able ji, where i is the argument passed to SplitDen. This routine is optimized 
for higher powers of the denominators and uses the sum formula rather than 
repeated splitting of pairs. 

SubsS 

Evaluates S functions with a numerical last argument into polynomials in the 
variables (if possible, with numerical values). For argument / = only nu- 
merical values are realized; for i = 1 the procedure also expands polynomials. 
S functions with symbolic last argument are untouched, of course. 

Sum It 

Calls the various procedures for summation algorithms. The argument speci- 
fies the index i of the sum to be done. 

SumPosPow 

Sums positive (or zero) powers of ji possibly in combination with one S or 
pow function. The argument specifies the index i of the sum to be done. 

SynchDenFac 

Synchronizes combinations of den and fac, invfac functions, i.e. factorials 
and denominators. The argument specifies the associated summation index /. 

SynchDenS 

Synchronizes combinations of den and S functions, i.e. denominators and 
S-sums. The argument specifies the associated summation index i. 

Synch FacS 

Synchronizes combinations of S and fac, invfac functions, i.e. factorials and 
5-sums. The argument specifies the associated summation index i. 

SynchNumPac 

Synchronizes combinations of (positive) powers of 7, and fac, invfac func- 
tions. The argument specifies the associated summation index i. 

SynchSS 

Synchronizes products of S'-sums. The argument specifies the associated sum- 
mation index /. 

As discussed in Section 2 the S-sums can be treated as generalization of the harmonic 
sums investigated in Refs. [1,9, 18]. It is therefore evident that some of the procedures 
presented here are similar to those in the Summer package [9] written by J. A.M. Ver- 
maseren. In particular the procedures AdjustSum, BasIsS, DoSynch, DoTheta, Flipj, 
Makej, SubsS, SynchDenFac and SynchSS are adapted versions of similar proce- 
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dures in SUMMER. 

3.2.3. Harmonic sums in infinity 

A special class of sums, occurring as a result of the summation algorithms, are har- 
monic sums in infinity. These are related to multiple zeta values [2], which for a given 
weight are reducible to a small set of transcendental numbers using, for example, the 
algebraic properties in Eq. (2.3). 

Together with the XSummer package, we provide the tables of harmonic sums in 
infinity (limited to weight six). These procedures are called tables, tablel , ... , tableG 

and the corresponding files are directly extracted from the SUMMER package [9]. This 
facilitates interfacing with routines of the current package. 

Let us also note that the reduction of multiple zeta values of a given weight to some 
irreducible set of constants is currently an active field of mathematical research. Cur- 
rently, downloads of tables up to weight nine in Form [23] or Maple format [24] are 
publicly available and extensions are known up to weight 16 [25, 26]. 

4. Examples 

Along with the distribution of the XSummer package, we provide also a number of 
non-trivial examples. These examples can either be run with the help of the shell 
script Testit in a standard Linux/Unix environment or with the help of TestltXP.bat 
under the Microsoft Windows XP operating system. (For either of the scripts the user 
might have to make small adaptations, though.) 

The respective script executes the Form files Examples.frm and Dolntegrals.frm (dis- 
cussed in detail below) and performs a check on the output of the former computations 
against tabulated results with the help of the file CheckResults.frm. In this way, the 
user can verify the correctness of the installation of the XSummer package. At the 

same time, the examples are meant to illustrate the usage of sums with the XSum- 
mer package. In particular, we want to clarify the conventions for the input to the 
procedures of Section 3. 

We provide in the file Examples.frm a number of (generalized) hypergeometric func- 
tions from the original Ref. [4], 

hypergeom2Fl (a*ep,b*ep,l-c*ep,xl) 
hypergeom2Fl (1, -ep, 1-ep, xl) 

hypergeom3F2 (-2*ep, -2*ep, 1-ep, 1-2 *ep, 1-2 *ep, xl) 
appel2 (1, 1, ep, 1+ep, 1-ep, xl, x2) 

where the coefficients of Laurent series in e are calculated up to a given order in terms 
of multiple poly logarithms. The first example hypergeom2Fl (a*ep, b*ep, l-c*ep, xl ) 
was actually given in Eq. (1.1). In Examples.frm it is realized by the following set of 
substitutions: 
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id hypergeom2Fl (a?,b?,c?,xl?) = suml ( jl, 0, inf ) * 

Po ( jl+a,a) *Po ( jl+b,b) *InvPo ( jl+c, c) *invfac( jl) * pow(xl, jl) ; 



id Po(xl?,x2?) = Gamma (xl) *InvGamma (x2) ; 
id InvPo (xl?,x2?) = Gamma (x2) *InvGamma (xl) ; 

Here Po (xl, x2) and InvPo (xl, x2) denote the Pochhammer symbols, see Eq. (1.1). 
Next, the file Dolntegrals.frm calculates certain one- and two-loop Feynman integrals 
used in a complete calculation of higher-order perturbative corrections in Quantum 
Chromodynamics [7]. To that end, integrals of various topologies had to be considered. 
The respective analytical representations in terms of nested sums valid for arbitrary 
powers of denominators for the so-called C-topology [4], the B-box (one-loop box with 
one external mass [27]) and one-loop triangles with up to two external masses [28] are 
all given in the procedure lnt2Sum. The input to Dolntegrals.frm, i.e. the information 
on the topology and on the particular values for the powers of denominators, is passed 
on by preprocessor variables in Form. The explicit lists of all inputs are contained in 
Selectlntegral. 

In the actual calculation, we choose dimensional regularization [29-32] with D = 4 — 
2e and the modified [33] minimal subtraction [34] scheme for all loop integrals in this 
section. For illustration we give here the explicit series representation as we use it in 
lnt2Sum for the one-loop triangle with two external masses Tri(m,Vi,V2,V3;xi). It is 
defined by 

Tri(m,Vi,V2,V3;^i) = (-.i2)-'"+^+"i^^ / ,2^vw , Ivs (^.1) 



and, in the following, we use the short-hand notation Vij = V; + Vy for sums of powers 
of propagators. We have k2 = k\—pi,k'i = k2— pi, and we define the quantities 

512=(P1+P2)^ x=^. (4.2) 

512 

Equation (4.1) can be written as a combination of hypergeometric functions 2^1. The 
series representation for this integral with |x| < 1 is given by 

r(e-m+V23)r(l-e + m-V23)r(m-£-Vi3) 
Tri(m,Vi,V2,V3;j:) = — 



i=0 ' 



r(vi)r(v2)r(v3)r(2m - 2e - V123) 

"^m-e-V23 r(/i+Vi)r(/i-e + m-V2) 
r(/i + l+m-8-V23) 

r(;i+V3)r(/i -/» + e + Vi23^ 
r(/i + l-m + e+V23) 



(4.3) 



Equation (4.3) is implemented in the procedure lnt2Sum and prepared for the input 
to the XSUMMER package with a set of substitutions similar to those briefly discussed 
above for the 2F1 hypergeometric function. 
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Now, with the examples at hand, we would like to discuss the efficiency of the imple- 
mentation of the XSuMMER package in Form. As already mentioned above, certainly 
one major disadvantage in using Form is the lack of internal algorithms for polynomial 
algebra. As a consequence, we cannot easily bring rational functions of polynomials 
to a normal form. For that purpose, we would need operations such as factorization, 
polynomial division or combinations thereof as used for instance in partial fractioning. 
In applications of the algorithms A-D to sums with multiple scales, this problem be- 
comes apparent when trying to normalize the arguments of the 5-sums. In particular, 
the recursions of the conjugation algorithm (type C) are sensitive to this issue. Here, 
efficient simplifications of polynomials can have a significant impact on the execution 
time for a given sum. 

Since no really elegant way for simplifying polynomials exists within Form, we have 
provided a case-by-case solution. Inside the procedure AlgC we do have the (user- 
accessible) procedure SimpleArgs (described in the previous section). There, depend- 
ing on the scales involved, appropriate substitutions for tuning or optimizations should 
be added. Then, we do have the procedure CustomizeDen. It uses partial fractioning 
and normalizes the ubiquitous denominators to standard factors. In a similar spirit, we 
also use the procedure ParFrac, which brings two-scale polynomials to normal form 
in Dolntegrals.frm by means of partial fractioning at the end of the calculation. 
Another more generic approach to this problem, which the experienced user might 
consider, is the following. FORM provides the preprocessor statements #system and 
#pipe, which invoke a call to the operating system. This provides the opportunity to 
realize a link to external programs, in particular to computer algebra systems such as 
Maple or Mathematica along with their functionalities in polynomial algebra. How- 
ever, both the detailed description and the implementation of these features is beyond 
the scope of the present work. 

Finally, we would like to give some information on the use of computer resources. 
Typically, the execution times and the use of memory or disk space depend very much 
on the problem under consideration. As a general rule, the expression size and there- 
fore the execution time correlate strongly with the depth of the Laurent expansion in 
8. In the examples, we can control this with the preprocessor variable EXPANDEP. This 
effectively cuts the series expansion in £ at the specified power. Note that for every 
individual term the series is cut to the specified order. If additional poles are present it 
might thus be necessary to expand individual terms to higher order to make sure that 
no terms are lost. For safeguard, we added a function order, which shows the effect 
of the truncation. Of course, the run time is also correlated with the number of nested 
sums. In practical applications it might be advantageous to start with a low value for 
EXPANDEP and increase it until the final result has the desired order. In addition to 
this, in practical applications, some tuning may be needed for new large problems, 
particularly in bringing polynomials to normal form. 

Let us close this section with a few remarks about the runtime performance. To give 
the interested user a hint on how long the presented examples will take, we measured 
the runtime on a standard PC. In particular, the machine we used was equipped with a 
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3 GHz P4 (hyperthreading support) with 512 MB of RAM. The considered problems fit 
completely into the RAM so that the hard-drive speed plays only a marginal role. The 
results are shown in Table 4, where we averaged always over 10 independent runs. 
The different Form versions used are available at http://www.nikhef.nl/~form. 







Dolntegrals.frm 


Form 3.1 executable 


Examples.frm 


CTOPO 


BBOX 


TRIMASS 


gcc, 15-jul-2005 


7.9s 


20.9s 


119.3s 


19s 


ice, 24-jan-2003 


7.5s 


19.9s 


108.2s 


17.4s 


gcc, 17-oct-02 


8.7s 


24.5s 


181.3s 


27.6s 



Table 2: Runtimes of the examples under Linux for different Form versions. 



We also tried a preliminary version of a native Microsoft Windows XP Form exe- 
cutable. The runtimes are very similar to those obtained with the icc-compiled FORM 
executable and the new gcc version. The difference of 5% might be due to the fact that 
on Windows XP the total time was measured and not only the user time. 
In applications more involved than the examples shown here, the user is advised to 
make a detailed analysis on the actual depth of the expansion in £ needed and on the 
particular structure of the polynomials. Furthermore the user should provide additional 
routines to simplify the symbolic arguments of S-sums. 

5. Conclusion 

Symbolic summation has advanced to an important method, e.g. in perturbative quan- 
tum field theory, and significant progress has been made during the past years. In the 
present work, we have provided an implementation in Form of algorithms suitable for 
the expansion of transcendental functions in a small parameter around integer values, 
such that the resulting (generalized) hypergeometric series can be expressed in closed 
form in multiple polylogarithms. 

As examples, we have discussed various applications in quantum field theory, particu- 
larly in the calculation of Feynman diagrams at higher orders in perturbation theory. In 
this context, the XSUMMER package has been used in perturbative calculations exten- 
sively and we believe it may be useful for a larger community. We have chosen Form 
for the implementation, because it is a fast and efficient computer algebra system and 
because of its capability to handle large expressions. For convenience of the user, we 
provide along with XSummer a set of sample calculations. These illustrate the use of 
the program. An extension of the present implementation to cover also algorithms for 
generalized sums from expansions around rational numbers (see e.g. [35]) will be the 
subject of a future publication. 
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Note added 

Very recently, Ref. [36] appeared, which addresses the problem of expanding hyper- 
geometric functions jFj^i around integer parameters to arbitrary order. It provides an 
implementation in Mathematica of the algorithms (2.7), (2.8), i.e. type A and B of the 
original Ref. [4]. Thus, it is capable of performing some expansions already discussed 
in Ref. [4] and discussed also in the present Letter, for instance Eq. (1.1). 
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